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ABSTRACT 

We present the results of a very deep (500 ks) Chandra observation, along with tailored nu¬ 
merical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc 
(19') to the northwest of M 87. The northern part of the front appears the sharpest, with a width 
smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere 
along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gra¬ 
dient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed 
across the interface. Such transport processes can be naturally suppressed by magnetic fields 
aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations 
indicates that in order to maintain the observed sharp density and temperature discontinuities, 
conduction must also be suppressed along the magnetic field lines. However, the northwestern 
part of the cold front is observed to have a nonzero width. While other explanations are pos¬ 
sible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) 
on length scales of a few kpc. Based on comparison with simulations, the presence of KHI 
would imply that the effective viscosity of the intra-cluster medium is suppressed by more 
than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent 
viscosity. Underneath the cold front, we observe quasi-linear features that are ~ 10 per cent 
brighter than the surrounding gas and are separated by ~ 15 kpc from each other in projection. 
Comparison to tailored numerical simulations suggests that the observed phenomena may be 
due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, 
where the magnetic pressure reaches ~ 5-10 per cent of the thermal pressure, reducing the 
gas density between the bright features. 

Key words: galaxies: X-rays: galaxies: clusters - galaxies: clusters: individual (Virgo) - 
clusters: intracluster medium - hydrodynamics - instabilities 


1 INTRODUCTION 

One of the first important results obtained with the Chandra X-ray 
Observatory was the surprising discovery of cold fronts - remark¬ 
ably sharp surface brightness discontinuities, with relatively dense 


and cool gas on the inner (bright) side and low density hot gas on 
the outer (faint) side. The first discovered cold fronts in Abell 2142 
and Abell 3667 were interpreted as signs of mergers ( |Markevitdi| 
|et al.|2000{|Vikhlinin et al.|2001a| l, where the sharp surface bright- 
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ness discontinuity separates the hotter intra-cluster medium (ICM) 
from the low entropy core of a merging subcluster. 

However, as more clusters were observed with Chandra, cold 
fronts were also discovered in the centers of cool-core clusters (at 
r =10-400 kpc) with no signs of recent major mergers ( |Markevitdi| 
|et al.|200r| |. Indeed, such cold fronts are observed in most relatively 
relaxed cool-core clusters ( [Ghizzardi et al.|2()T0] |, often with sev¬ 
eral, oppositely placed arcs curved around the central gas density 


peak, forming a spiral-like pattern (see e.g. Markevitch & Vikhlinin 


2007 

[Owers et al. 

[2009 [Simionescu et al.[[2012| [Rossetti et al. 

2013 

[Ichinohe et a] 

.|2015|. Considering projection effects and the 


weak density jumps, possibly all cool-core clusters may host some 
cold fronts. 

Cold fronts in cooling cores are believed to be due to the slosh¬ 
ing (or swirling) of the gas in the gravitational potential of the clus¬ 
ter. Numerical simulations show that sloshing can be induced easily 
by a minor merger where a subcluster falls in with a nonzero im¬ 
pact parameter ( [Tittley & Henriksen|2005[[Ascasibar & Markevitch| 
|2006| ). Due to this disturbance, the density peak of the main cluster 
swings on a spiral-like trajectory relative to its center of mass. As 
the central ICM is displaced from the corresponding dark matter 
peak, opposite and staggered cold fronts are created where cooler 
and denser parcels of gas from the centre come into contact with 
the hotter gas at larger radii. The large-scale coherent flows associ¬ 
ated with sloshing and the corresponding cold fronts can persist for 
Gyrs. By redistributing the ICM, sloshing plays an important role 
in shaping the cluster cores and may even help prevent significant 
cooling fiows ( [ZuHone et al.|2010]|Keshet|2012| l. 

|Markevitch et al. | ( [200 1 1 ) showed that sloshing cold fronts re¬ 
quire significant acceleration of the gas beneath the contact dis¬ 
continuity. This can be explained as centripetal acceleration of a 
tangential bulk flow just below the cold front ( [Keshet et al.|2010| l, 
which causes a strong tangential shear along the discontinuity. This 
strong shear flow should promote the growth of hydrodynamic in¬ 
stabilities, which would broaden and eventually destroy the cold 
front on short sub-dynamical time scales. But cold fronts are re¬ 
markably sharp, both in terms of the density and the temperature 
jumps: the gas density discontinuity in e.g. the merging cluster 
Abell 3667 is at least several times narrower than the Coulomb 
mean free path ( [Vikhlinin et al.|[2001a[ [Markevitch & Vikhlininj 
[2007[ ) and the observed temperature jumps require thermal conduc¬ 
tion across cold fronts to be suppressed by a factor of order ~ 10^ 
compared to the collisional Spitzer or saturated values ([Ettori &[ 
[Fabian|2000|[Xiang et al.|2007] t. A layer of magnetic field parallel 
to the discontinuity could stabilize the front, prevent the growth of 
instabilities, and suppress transport across the cold front ([Vikhlinin[ 
[et al.|2001b[ l. 

[Vikhlinin et al.[ ( [2001b[ ) and [Vikhlinin & Markevitch] ( [2002[ ) 
used the stability of the merger induced cold front in Abell 3667 
to derive a lower limit on the magnetic field of B >7-16 jiG at 
the discontinuity. Magnetic amplification in sloshing cold fronts is 
achieved through the shear due to the bulk flow of plasma beneath 
the front, which naturally sustains the strong magnetic fields par¬ 
allel to cold fronts ( [Keshet et aLpOlO) [ZuHone et aLpOlI) . In 
the layers underneath cold fronts, fields with initial strengths of 
0.1-1 jdG may be amplified to tens of jdG, and the magnetic pres¬ 
sures may be amplified to values reaching a significant fraction of 
the thermal pressure, leading to gas depleted bands ([ZuHone et al.[ 
1 ^ . 

One of the most interesting aspect of cold fronts is that they 
provide a relatively clean experimental setup. This allows us to 
address questions about the physics of the ICM in a more exact 



Figure 1. XMM-Newton EPIC/MOS mosaic image of the central r ~ 
150 kpc region of the Virgo cluster. The exposure corrected image has been 
divided by the best-fitting radially symmetric beta model (see [Simionescu[ 
[et al.|2010j[Roediger et al.|2()TT} . The white square indicates our Chandra 
ACIS-I pointing shown in Fig. 2. The “+” sign indicates the center of M 87 
at the bottom of the global gravitational potential of the Virgo Cluster. The 
“X” sign indicates the center of curvature (RA: 187.676, DEC: 12.5064) of 
the part of the cold front imaged by our Chandra pointing. 


way than in cluster cores with complex astrophysical substructure. 
The closest example of a cluster cold front is the spectacular sur¬ 
face brightness discontinuity that lies 90 kpc (19') to the north of 
M 87 in the Virgo Cluster and extends 140° in azimuth (see Fig.[^ 
Simionescu et al.|2010|). At the distance of 16.1 Mpc ([Tonry et 5![ 
2001> , 1 arcsec corresponds to only 78 pc, making this the best re¬ 
solved cluster cold front in the sky. The gas density jump across 
the discontinuity is Pm/Pout ~ 1.4 and the gas in the inner, bright 
side is ~40% more abundant in Fe than the ICM outside the front 
( [Simionescu et al.|2010| l. [Simionescu et al.[ ( [2007[ t also detected a 
much weaker cold front ~ 17 kpc to the south-east of the core. The 
opposite and staggered placement of these fronts, as well as the ab¬ 
sence of a remnant sub-cluster, make gas sloshing the most viable 
explanation for the presence of these edges ( [Tittley & Henriks^ 
[2005|[Ascasibar & Markevitch[2006[ ). 

Using constrained hydrodynamical simulations, [Roediger| 
[etal.[ ( [20TT| l showed that a sloshing scenario can reproduce the radii 
and the contrasts in X-ray brightness, projected temperature, and 
metallicity across the observed cold fronts in the Virgo Cluster. By 
comparing synthetic and real observations, they estimate that the 
original minor merger event that triggered the sloshing took place 
about 1.5 Gyr ago, when a sub-cluster of 2-4 x 10^^ M© passed the 
Virgo core at 100 to 400 kpc distance nearly perpendicular to our 
line-of-sight. The resulting gas sloshing, with motions primarily in 
the plane of the sky, creates a contact discontinuity with a surface 
that is predominantly perpendicular to our line of sight, maximizing 
the contrast and minimizing projection effects. In subsequent work, 
[Roediger et al.[ ( [20T^ demonstrated the capability of deep Chandra 
observations to reveal Kelvin-Helmholtz instabilities in the Virgo 
cold front by detecting multiple edges and shallower than expected 
surface brightness profiles across the discontinuity. Based on rela- 
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Table 1. Summary of the Chandra observations. Columns list the observa¬ 
tion ID, observation date, and the exposure after cleaning. 


Obs. ID 

Obs. date 

Exposure (ks) 

15178 

2014 Feb. 17 

47.1 

15179 

2014 Feb. 24 

41.9 

15180 

2013 Aug. 1 

140.5 

16585 

2014 Feb. 19 

45.6 

16586 

2014 Feb. 20 

49.8 

16587 

2014 Feb. 22 

37.8 

16590 

2014 Feb. 27 

38.1 

16591 

2014 Feb. 27 

23.8 

16592 

2014 Mar. 1 

36.1 

16593 

2014 Mar. 2 

38.1 


lively short XMM-Newton observations, they also present tentative 
indications that the northeastern portion of the front is disturbed by 
instabilities. 

Here, we present a very deep (500 ks) Chandra observation of 
the sharpest, northwestern part of the Virgo cold front, designed to 
test the theoretical predictions in unprecedented detail. 

2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Chandra data 

The Chandra X-ray observatory observed the prominent cold front 
to the northwest of M 87 (see Fig. between 2013 August 
and 2014 March using the Advanced CCD Imaging Spectrome¬ 
ter (ACIS-I). The standard level-1 event lists were reprocessed us¬ 
ing the CIAO (version 4.7, CALDB 4.6.5) software package, in¬ 
cluding the latest gain maps and calibration products. Bad pixels 
were removed and standard grade selections applied. The data were 
cleaned to remove periods of anomalously high background. The 
observations, along with the dates, identifiers, and net exposure 
times after cleaning (total net exposure of 500 ks), are listed in 
Table 1. Background images and spectra were extracted from the 
blank-sky fields available from the Chandra X-ray Center. These 
were cleaned in an identical way to the source observations, repro¬ 
jected to the same coordinate system and normalized by the ratio 
of the observed to blank-sky count rates in the 9.5-12 keV band. 

2.1.1 Image analysis 

Background-subtracted images were created in nine narrow en¬ 
ergy bands (0.8-L0 keV, 1.0-1.2 keV, 1.2-1.5 keV, L5-2.0 keV, 
2.0-2.5 keV, 2.5-2.75 keV, 2.75-3.0 keV, 3.0-3.5 keV, and 3.5- 
4.0 keV), spanning 0.8-4.0 keV. The nine narrow-band images 
were fiat fielded with respect to the median energy for each image 
and then co-added to create the X-ray images shown in Fig.|^ In the 
0.8-4.0 keV band the signal over background ratio is optimal and 
the emissivity is nearly independent of the gas temperature, max¬ 
imizing the contrast at which we detect the contact discontinuity. 
Identification of point sources was performed using the CIAO task 
WAVDETECT. The point sources were excluded from the subsequent 
analysis. 

In order to emphasize structure in the image, we also produced 
a beta-model subtracted image, shown in the right panel of Fig. 

In order to remove the underlying large-scale surface brightness 
gradient, we fitted a beta-model to an azimuthally averaged X-ray 
surface brightness profile extracted from the XMM-Newton mosaic 
image (see Sect. |2.2| and Fig. and divided our Chandra image 


by this model, patched on large scales with an 80 arcsec smoothing 
window (see |Zhuravleva et al.|2015| for details). Our patched model 
with which we divide our image is defined as /pm = IpSa-Ux! 
where Ip is the best fit beta-model to the azimuthally averaged X- 
ray surface brightness profile, S^rL] is Gaussian smoothing with 
window size cr, and lx is the X-ray surface brightness image. 

2.1.2 Spectral analysis 

Individual regions for the 2D spectral mapping were determined 
using the contour binning algorithm ( |Sanders|2006| l, which groups 
neighbouring pixels of similar surface brightness until a desired 
signal-to-noise ratio threshold is met. We adopted a signal-to-noise 
ratio of 75 (approximately 5600 counts per region), which provides 
approximately 3 per cent precision for density, 4 per cent precision 
for temperature, and 10-20 per cent precision for metallicity mea¬ 
surements, respectively. We modelled the spectra extracted from 
each spatial region with the SPEX package ( [Kaastra et al.||1996] l. 
The spectral fitting was performed in the 0.6-7.5 keV band. The 
spectrum for each region was fitted with a model consisting of an 
absorbed single-phase plasma in collisional ionization equilibrium, 
with the spectral normalization (emission measure), temperature, 
and metallicity as free parameters. The deprojected spectra are ob¬ 
tained following the method described by |Churazov et al.|{|2003|) 
and are fitted using XSPEC and the APEC plasma model ( [Smith et al.| 
|2001 [[Foster et al.|2012| |. The line-of-sight absorption column den¬ 
sity was fixed to the value. Ah = 2 x 10^^ cm“^, determined by the 
Leiden/Argentine/Bonn radio survey of Hi ( [Kalberla et al.|2003] ). 
Metallicities are given with respect to the proto-Solar abundances 
by [Lodders[ ( |2003l ). 

2.2 XMM-Newton data 

A mosaic of six XMM-Newton pointings, each with an exposure 
time of ~ 20 ks, was performed in June 2008. An additional point¬ 
ing was obtained in November 2008. These observations cover an 
annulus around the deep (net exposure time of 120 ks) pointings 
on M 87. The full exposure corrected mosaic of XMM-Newton ob¬ 
servations divided by the best fit radially symmetric beta-model is 
shown in Fig.[^ See [Simionescu et n^ ( [2010[ ) for the detailed dis¬ 
cussion on the data and the analysis. 


3 RESULTS 

The most prominent sloshing cold front in the Virgo Cluster ex¬ 
tends at r ~ 90 kpc north of M 87 (see Fig.[^. Our Chandra ob¬ 
servation (indicated by the white square in the figure) targets the 
part of the front with the largest surface brightness contrast. The 
deep Chandra image, on the left panel of Fig.|^ shows clearly this 
surface brightness discontinuity. 

The Chandra image also reveals unexpected, relatively bright 
quasi-linear features underneath the discontinuity. To our knowl¬ 
edge, this is the first time that such substructure has been identified 
with a cold front. Three of these quasi-linear features can be seen 
clearly on the residual image on the right panel of Fig.[^ They are 
~ 10 per cent brighter than the surrounding gas and are separated 
by ~ 15 kpc from each other in projection. These features are below 
the detection limit of the existing XMM-Newton observation. 

To study the detailed cold front interface, we extracted sur¬ 
face brightness profiles along 15 narrow, 5 degree wedges, each 
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Figure 2. Left panel: Chandra image of the cold front in the 0.8-4.0 keV energy band. The image was smoothed with a Gaussian function with a 1.5 arcsec 
window. Right panel: The residual image obtained by dividing the image on the left with a beta-model, and patched on large scales with an 80 arcsec scale 
smoothing window (see |Zhuravleva et al.|2015| and the text for details), reveals three X-ray bright quasi-linear features separated from each other by ~ 15 kpc 
in projection. The dark band outside the cold front is an image processing artifact due to the sharp surface brightness discontinuity. Point sources were removed 
from the residual image. 


corresponding to ~ 5 kpc along the front (see Fig. [^. To opti¬ 
mize the sharpness of the surface brightness edge in our profiles, 
we aligned our annuli with the front interface, choosing the ver¬ 
tex of the wedges at the center of curvature of the highest-contrast 
part of the cold front (RA: 187°.676, DEC: 12°.5064). The data 
were binned into 4.92 arcsec radial bins. We fit the surface bright¬ 
ness profiles by projecting a spherically symmetric discontinuous, 
broken power-law density distribution (see e.g. |Owers et al.|2009 


The fitting was performed using the PROFFIT package ( Eckert et^ 
|2011| ) modified by |Ogrean et al.| ( [T015| ). We fit the data in the ra¬ 
dial range of 8-16 arcmin, corresponding to 38-75 kpc from the 
center of curvature. The free parameters of our model are the nor¬ 
malization, the inner and outer density slopes, and the amplitude 
and radius of the density jump. We verified that the best fit parame¬ 
ters are not sensitive to the utilized radial range. The surface bright¬ 
ness profiles over-plotted with their best fit models, along with the 
values for the best fit density jumps, radii of the discontinuity, and 
reduced;^f^s (for 92 degrees of freedom) are shown in Fig.[^ 


this resolved gradient, our upper limit on the width is significantly 
larger than Chandra’s resolution. 

On the other hand, at azimuths of 25-70 degrees, our model 
mostly provides a poorer fit to the data. Fig.j^shows that around the 
brightness edge most profiles display deviations from the model, 
with shallow, smeared discontinuities. The parameters of the best 
fit model and the structure seen in the profiles also show azimuthal 
variation. In particular, the northwestern part of the front, at az¬ 
imuths below 55 degrees, shows significant variation both in the 
value of the best fit density jump and radius. For the sectors with 
the least sharp edges the fit significantly improves after including 
Gaussian smoothing in the model. For example, for the profile ex¬ 
tracted in the 65-70 degree wedge, Gaussian smoothing with a best 
fit width of 4.3 kpc improves the from 138 (for 92 degrees 
of freedom) to 101, for one additional free parameter. According 
to the F-test, the likelihood that this is a chance improvement is 
1.06 X 10"^. This highly significant improvement strongly suggest 
that at these azimuths the front has a finite width. 


Our model provides the best fit to the northern part of the front, 
at azimuthal angles of 70-95 degrees, indicating that this is the 
sharpest part of the cold front. In particular, our model provides the 
best fit to the profile extracted at the azimuthal range of 80-85 de¬ 
grees. In order to place an upper limit on the width of the front, we 
fit this profile with a model that includes Gaussian smoothing as an 
additional free parameter. We find that the 99 per cent upper limit 
on the smearing of this part of the cold front is 2.5 kpc, which is 
1.5 Coulomb mean free paths (see Sect. |5.1| for details). Given the 
quasi-spherical geometry along our line of sight, in this nearby sys¬ 
tem we observe a steep, few kpc wide surface brightness gradient at 
the front (in the more distant systems where the gradient cannot be 
resolved we only see a sharp surface brightness jump). Because of 


To study the projected thermodynamic properties and the 
metallicity of the ICM across the cold front, we extracted spectra 
along three 22.5 degree wide wedges (labeled North, Northwest, 
and West), the vertex of which is at the same center of curvature of 
the cold front as for the surface brightness profiles. The best fit pro¬ 
jected temperature and metallicity profiles along the three wedges, 
determined in narrow 2 kpc wide regions, are shown in Fig. 
The projected temperature as a function of radius increases from 
~ 2.7 keV to ~ 3.6 keV across the interface, and the metallicity 
decreases from ~ 0.6 Solar to ~ 0.3 Solar. The projected temper¬ 
ature distribution appears to be smeared across the cold front. The 
observed gradient of the projected temperature is between ~ 4 kpc 
(West) and 8 kpc (North) wide. 
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Figure 3. Surface brightness profiles across the cold front measured along 15 narrow, 5 degree wedges, each corresponding to ~ 5 kpc along the front. The 
azimuth angles are measured counterclockwise from the west. Radii are measured from the center of curvature of the cold front (see Fig.[^. Over-plotted are 
the best fit broken power-law models for the ICM density profiles. For each wedge we indicate the best-fit radius of the discontinuity (Rf, measured from the 
center of curvature of the cold front), density jump (nj), and the reducedfor 92 degrees of freedom. 


The 2D maps of the projected thermodynamic properties 
shown in Fig. [^confirm these trends. While the temperature clearly 
increases across the front, the best fit metallicity drops from ~ 0.6 
Solar to ~ 0.3 Solar. The projected temperature and entropy distri¬ 
butions appear slightly smeared at the edge. The smearing is likely 
to be at least in part due to projection effects (see Sect. |5.1| for more 
discussion). As expected, the cold front is not visible in the pres¬ 
sure map, indicating that the sloshing gas is in approximate hydro¬ 
static equilibrium within the cluster potential. The linear structures 
identified in the images (see Fig.[^ do not show any significant dif¬ 


ference in temperature or metallicity from the surrounding ambient 
medium. 

In order to study the temperature difference of the quasi-linear 
features and the ambient plasma in more detail, we extracted five 
spectra from rectangular regions: three spectra from regions cen¬ 
tered on the quasi-linear features, and two spectra from the sur¬ 
rounding regions. Based on fitting the combined spectra, the best 
fit temperature of the quasi-linear features is kT = 2.77 + 0.03 keV 
and the best fit temperature of the surrounding gas is kT = 2.92 + 
0.03 keV, qualitatively consistent with the trend expected from the 
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Figure 4. Projected temperature and metallicity profiles along three 22.5 
degree wide wedges, extracted from the northern direction westward. The 
profiles extracted from the north, northwest, west are shown in black, blue, 
and red, respectively. The best fit radius of the front determined by fitting the 
surface brightness profiles is indicated by the dashed line (the uncertainty on 
the front radius is +2 kpc). Radii are measured from the center of curvature 
of the cold front. 
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simulations we present in Sect.|^(see Fig.|^. However, we caution 
that because the surface brightness excess associated with these 
features is only ~ 10 per cent, the measured temperature differ¬ 
ence could easily be due to the underlying spatial variation in the 
temperature of the ambient gas. 

In order to determine the deprojected spectral properties of 
the ICM across the cold front, we also extracted spectra from a 
wedge centered on M 87 at the bottom of the gravitational poten¬ 
tial of the Virgo Cluster. Centering the series of partial annuli used 
for spectral extraction at the bottom of the global gravitational po¬ 
tential enables us to perform the most robust deprojection of the 
stratified cluster atmosphere. The 30 degrees wide wedge used for 
extracting the Chandra spectra extends from the north towards the 
west. For this particular wedge, the width of the cold front is the 
least affected by the choice of the different center. Fig. shows 
the comparison of the deprojected electron density, temperature, 
metallicity, pressure, and entropy profiles determined in the wedge 
using Chandra data (red datapoints) with the properties determined 
in the other directions (averaged over a 330 degree region that ex¬ 
cludes the wedge) using XMM-Newton spectra (blue datapoints). 
While the deprojected pressure appears continuous across the cold 
front, the density, temperature, entropy, and metallicity distribu¬ 
tions show clear discontinuities. The deprojected entropy profile 
extracted across the cold front shows a remarkably flat distribution 
underneath the discontinuity that extends across a radial range of 
~ 50 kpc. The comparison also shows that under the cold front 
(~ 20 kpc from the discontinuity) the ICM entropy is smaller and 
above the front it is higher than in the other directions. The pressure 
appears slightly higher under the cold-front than in the other direc¬ 
tions. We verified that the deprojected profiles determined using 


Figure 5. Deprojected electron density, temperature, metallicity, pressure, 
and entropy profiles (from the bottom to the top), respectively. The red dat¬ 
apoints were determined in a 30 degree wedge extending from the north to 
the west, centered on M 87, using Chandra data and the blue datapoints 
are averaged over the other directions (330 degree region that excludes the 
wedge) using XMM-Newton data. The radius of the cold front is indicated 
with the dotted line. While the deprojected pressure is continuous across 
the cold front, the density, temperature, entropy, and metallicity distribu¬ 
tions show relatively strong discontinuities. Note that under the cold front 
(~ 20 kpc) the ICM entropy is smaller and above the front it is higher than 
in the other directions. 

Chandra and XMM-Newton data in the 30 degree wedge show con¬ 
sistent results. We caution that while our deprojection procedure 
assumes spherical symmetry centered on M87, the distribution of 
the gas uplifted by sloshing is likely to be different and therefore 
our inferred thermodynamic properties underneath the discontinu¬ 
ity may be biased. 


4 NUMERICAL SIMULATIONS 

We compare our long exposure of the Virgo cold front to mag¬ 
netohydrodynamic (MHD) simulations. These simulations were 
performed using the Athena 4.1 astrophysical simulation code 
( [Stone et al.|2008|). The tech nical details of the simulation setup are 
discussed in jZuHone et al.| ( |20I5| ) and the initial conditions, based 
on the setup from jRoediger et al^ ( 201 Ij ) and jRoediger et al.| ( |2013| ), 
were chosen to match the global Virgo Cluster potential and the 
ICM properties. The cluster is modeled as a cooling-core system 
with a mass of ~ 2 X 10^"^ M© and ~ 2 keV. The sloshing in the 
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Figure 6. Projected thermodynamic maps. Each spatial region contains 5600 net counts (S/N of 75). To compare the pseudo-pressure and pseudo-entropy dis¬ 
tributions, we assume the gas uniformly distributed along the line-of-sight depth of / = 1 Mpc over the entire field of view. The units of temperature (top left), Fe 
abundance (top right), pressure (bottom left) and entropy (bottom right) are keV, Solar jLodders|2003) , 1000 keV cm“^(//l Mpc)^/^ and keV cm^(//l Mpc)^/^, 
respectively. The fractional Icr statistical errors are ~ 4 per cent for temperature, 5 per cent for pressure and entropy, and ~ 10-20 per cent for metallicity. 
Point sources were excluded from the regions shown in grey. 


simulation is due to the passage of a gas-less dark-matter subcluster 
with a mass of 2x10^^ M© at r ~ 100 kpc at the time of r ~ 1 Gyr af¬ 
ter the beginning of the simulation. We note that simulations show 
that similar cold front features are produced by gas-less perturbers 
approaching at close distances or gas-filled perturbers approaching 
at larger distances, indicating that the choice of perturber is not as 
important as the shape of the potential well and the thermodynamic 


profiles of the main cluste r {[ZuHone et aL||2010[ |Roediger et al.| 
|2011| |Roediger & Zuhone||2012| l. The time at which the proper¬ 

ties of the Virgo cold front most resemble the results of the sim¬ 
ulation is ~ 1.7 Gyr after the passage of the subcluster. In Fig.[^ 
we show the synthetic X-ray surface brightness residual distribu¬ 
tions for the simulation at this time for two pure hydrodynamic runs 
with different viscosities (top row), and for two MHD runs (bottom 
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row). The pure hydrodynamic simulations were performed assum¬ 
ing either inviscid ICM {ji - 0), or 10 per cent of the isotropic 
Spitzer-like temperature dependent viscosity (// = O.lyUsp). Both 
MHD runs started with a tangled magnetic field distribution (set up 
as in [ZuHone et al.|201l] >, one with an initial ratio of thermal to 
magnetic pressure of p - 1000 and the other with p - 100. 

The two pure hydrodynamic simulations are visibly inconsis¬ 
tent with the data: the inviscid run shows large, well-developed 
Kelvin-Helmholtz instabilities (KHI) that are not seen at the Virgo 
cold front image (the surface brightness profiles extracted from the 
mock image show a distorted front with multiple strong edges; see 
also Fig. 4 in |Roediger et al.|2013] t; and the run with 10 per cent 
of Spitzer viscosity produces a discontinuity that is much smoother 
than the profiles extracted across the northwestern part of the cold 
front, seen in Fig. Importantly, pure hydrodynamic simulations 
do not reproduce the bright linear features seen in the right panel of 
Fig. H though such features appear in the MHD run with the initial 
plasma p = 100; we will discuss this in Section [53] 


5 DISCUSSION 


taking place predominantly in the plane of the sky projection ef¬ 
fects are minimized. Their smearing effect gets more important as 
the shape of the front departs from spherical symmetry and as the 
inclination of the orbital plane to the plane of the sky gets larger 
(until the inclination approaches 90 degrees, when sloshing is along 
the line-of-sight). 

Unsuppressed diffusion would smear the density discontinuity 
by several mean free paths. Our upper limit on its intrinsic width 
therefore indicates that diffusion is suppressed across the cold front. 
This is consistent with the results obtained for the cold front identi¬ 
fied with a subcluster merger in Abell 3667, where [^khlinin et al.j 
( [2001 a| ) showed that the gas density discontinuity is several times 
narrower than the Coulomb mean free path (see [Markevitch &| 
|Vikhlinin|2007| for updated results). [Churazov & InogamoT { |2004j l 
argue that a small but finite width of the interface set by diffusion, 
that is of the order of a few per cent of the curvature radius (~ 2 kpc 
for the Virgo cold front), could strongly limit the growth of pertur¬ 
bations such as KHI (see also |Chandrasekhar|I96I| l. Given projec¬ 
tion effects that are expected to be broadening the observed width 
of the discontinuity, the intrinsic width of the Virgo cold front is 
likely to be too small to be suppressing the development of KHI. 


5.1 The width of the cold front and suppressed diffusion 


The observed surface brightness profiles across the cold front show 
considerable azimuthal variation. While at certain azimuth angles 
(particularly in the north) the profiles are relatively well described 
by the model of a spherically symmetric cloud with a power- 
law density distribution, in other directions they deviate from the 
model, are wider, and the best fit front parameters show significant 
azimuthal variation. 

For the best fit deprojected temperatures and electron densi¬ 
ties, the Coulomb mean free paths inside and outside of the cold 
front are Tin = 0.77 kpc and Tout = 2.2 kpc, respectively. For un¬ 
suppressed diffusion, particles diffusing from the bright, dense side 
of the discontinuity to the low density outer ICM will be especially 
efficient at broadening the cold front. The Coulomb mean free path 
for diffusion from the inner bright side to the outer faint side of the 
discontinuity is: 


din- 


>out 


( Tout \ / ^e,out 

7 keV / \ 10"^ cm“^ 


Tin \ G(l) 

Tout/G(V7W^) 


1.7 kpc. 


( 1 ) 

where We,out is the electron density outside the cold front, Tin and 
Tout are the temperatures inside and outside the cold fr ont, and 
G{x) = [(l)(x)-x(l)'{x)]/2x^, where (p(x) is the error function ( jSpitzer 
1 1962 [[Markevitch & Vikhlinin|2007| l. Our conservative upper limit 
of 2.5 kpc on the intrinsic width of the cold front, determined in the 
sector with the sharpest brightness edge, corresponds to about 1.5 
Coulomb mean free paths. 

Part of the apparent broadening of the cold front in the Virgo 
Cluster may be caused by projection effects and the true intrin¬ 
sic width of the discontinuity is also likely to be smaller than the 
Coulomb mean free path. Because the front is seen along the sur¬ 
face in projection, any deviations from the ideal spherical shape 
will smear the edge. Therefore, the observed width is a conserva¬ 
tive upper limit on the actual width of the discontinuity. Based on 
comparisons to numerical simulations, [Roediger et al.[ ( [201 1[[201^ 
argue for sloshing in the Virgo Cluster close to the plane of the 
sky, ruling out orbital planes that are further than 45 degrees from 
the plane of the sky. We verified that the mock X-ray images of 
sloshing in our simulations, produced for various lines of sight, are 
consistent with those in [Roediger et al.[ (2011 [[201^ . For sloshing 


5.2 Kelvin-Helmholtz instabilities and ICM viscosity 

The primary effect of projection is the broadening of the surface 
brightness edge and possibly a relatively smooth azimuthal change 
of the broadening. However, the non-zero width of the cold front, 
combined with the observed substructure and the azimuthal varia¬ 
tions of the best fit radius and density jump seen in Fig. [^ require 
the presence of additional processes, such as KHI. The temperature 
profiles shown in Fig. [^ are consistent with the presence of insta¬ 
bilities on a scale of a few kpc. 

Constrained hydrodynamic simulations by [Roediger et al.| 
( [2012[[20T3] | demonstrate the efficiency of viscosity at suppressing 
KHI and show that for viscosities of about 10 per cent of the Spitzer 
value or larger, at this cold front, KHIs should be suppressed. Our 
pure hydrodynamic simulation run with 10 per cent Spitzer viscos¬ 
ity is consistent with these conclusions. Based on comparison to 
these numerical simulations, the presence of KHI would imply that 
the effective viscosity of the ICM is suppressed by more than an 
order of magnitude with respect to the isotropic Spitzer-like tem¬ 
perature dependent viscosity. Our conclusions, however, depend on 
the correctness of the assumed merger model. 

The growth of KHIs may be facilitated by the presence of 
magnetic fields which lower the effective viscosity of the plasma 
by preventing the diffusion of momentum and heat perpendicular to 
the field lines. Also, micro-scale MHD plasma instabilities may set 
an upper limit on the viscosity that is much lower than that expected 
from ion collisionality ( [Kunz et al.|2014[ ). However, depending on 
the field geometry and strength, magnetic fields may also protect 
cold fronts against instabilities. A layer of strong enough magnetic 
fields aligned with the cold front could suppress the growth of KHIs 
( [Vikhlinin & Markevitch|200^ [Keshet et al.|2010[ ). However, the 
suppression of KHI in the simulations of [ZuHone et al.[ ( [2011[ ) is 
often partial, indicating that magnetic fields cannot always stabi¬ 
lize cold fronts. Our example in the bottom-right panel of Figure [^ 
shows that if magnetic fields are initially weak (p = 1000 in our 
example), they will be unable to suppress KHI. Stronger magnetic 
fields (see bottom-left panel of Fig. |3for an example with initial 
P - 100) suppresses most of the KHI but not all. 

Other possible explanations for the broadening and the ob¬ 
served substructure in the surface brightness profiles include 
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Figure 7. Simulated X-ray surface brightness residual distributions for the Virgo Cluster simulation att = 1.7 Gyr after the closest passage of the subcluster 
for two pure hydrodynamic simulation runs with different viscosities (top row), and two MHD runs with different initial plasma beta parameters (bottom row). 




Figure 8. Mock surface brightness images (left two panels), and surface brightness and temperature profiles (right two panels) extracted from the black wedges 
for numerical MHD simulations performed assuming no conduction and full Spitzer conduction along the magnetic field lines. Spitzer conduction along the 
magnetic field lines would wash out the temperature discontinuity to a degree that is inconsistent with the observations. The observed temperature and surface 
brightness profiles (see Fig.|^[^ are consistent with the conduction along the magnetic field lines being suppressed by a factor >10. 
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10 Wemer et al. 



Figure 9. A zoom in on a 80 x 80 kpc^ region around the cold front located northwest of the cluster center for the MHD simulation run with an initial magnetic 
to thermal pressure ratio = 100 (see the bottom right panel of Fig.^. Top-left: slice of gas density; top-right: slice of gas temperature; bottom-left: slice 
of magnetic field strength; bottom-right: projected surface brightness residuals. The solid cyan lines indicate the extraction regions for the profiles shown in 

Fig.Iig 


plasma depletion due to regions of amplified magnetic fields be¬ 
low the cold front, jet inflated bubbles entrained by the sloshing 
gas, or clumpy, inhomogeneous ICM distribution. 

5.3 Gas mixing at the cold front 

Well developed, long lasting, large KHI would mix the high metal- 
licity ICM (~ 0.6 Solar) on the bright side of the interface with the 
low metallicity gas (~ 0.3 Solar) on the faint side. The relatively 
sharp (narrower than 6 kpc) metallicity gradient indicates that if 
the surface brightness substructure is due to KHI, they are on rela¬ 
tively small, few kpc, scales and are unable to efficiently mix met¬ 
als across wider regions around the interface. Since cold fronts are 
transient wave phenomena, the lack of significant mixing across the 
front is not entirely surprising. Numerical simulations indicate that 
metallicity jumps remain preserved even for inviscid ICM ( |Roedi-| 
|ger et al.|2011| l and they get washed out only on very small scales, 
consistent with our observation and observations of sloshing cold 


fronts in other more distant clusters which are also associated with 
metallicity discontinuities (see e.g. |Ghizzardi et al.|2014| |. 

5.4 Conduction across the discontinuity 

|ZuHone et al.| < |20I3] l show that while sloshing approximately aligns 
the magnetic field lines with the cold front surface, perfect align¬ 
ment is not expected. KHI will re-tangle the field lines, restoring the 
heat flow between the gas phases above and below the cold front. 
|ZuHone et al.| ( |2015| ) also performed simulations for the Virgo 
Cluster with anisotropic conduction, assuming full Spitzer conduc¬ 
tivity along the magnetic field lines. In these, the temperature and 
surface brightness gradients are washed out to such a degree that 
is inconsistent with our observation (see Fig|^. The relatively nar¬ 
row temperature gradients across the interface (see Fig.|^ suggest 
that the effective conductivity along the field lines must be reduced 
by effects not modeled by the simulation, such as strong curva¬ 
ture of the field lines, stochastic magnetic mirrors on small scales. 
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Figure 10. Profiles of electron density, magnetic field strength, and pressure along the cyan lines indicated in Fig.[^as ‘1’ (left panel) and ‘2’ (right panel). The 
difference between the total and thermal pressure is due to magnetic pressure support. Dashed lines indicate approximate boundaries of density enhancements 
along the profile. 


and/or microscopic plasma instabilities (e.g. Chandran et al.|199'^ 


Malyshkin|2001|[Narayan & Medvedev|2001[|Schekochihin et al.| 

2005||2008| ). While the observed temperature and surface bright¬ 
ness profiles are consistent with no conduction along the magnetic 
field lines (see Fig.[^, conductivity at or below 0.1 of the Spitzer 
value along the field lines (which would still result in a steep gra¬ 
dient across the discontinuity, see |ZuHone et al.|2013| ) cannot be 
excluded and the exact suppression factor remains unknown. Like 
our conclusions on viscosity, this result also depends strongly on 
the correctness of the assumed merger model. 


5.5 Amplified magnetic fields underneath the cold front 

The Chandra images in Fig. show quasi-linear features that are 
~ 10 per cent brighter than the surrounding gas. The simulations of 
[ZuHone et al.| ( |20lf] ) and |ZuHone et al.|f2015| ) show that the mag¬ 
netic field layers amplified by sloshing may be wide and extend rel¬ 
atively far below cold fronts. This can be seen clearly in the Virgo 
simulation with the initial plasma p - 100 in Fig. [9] which shows 
slices through a 80 X 80 kpc^ region around the cold front for gas 
density, temperature, and magnetic field strength, as well as resid¬ 
uals of projected surface brightness. Along the lines indicated as 
“1” and “2”, we calculate the distributions of the density, tempera¬ 
ture, magnetic field strength and pressure and show them in Fig.p^ 
The quasi-linear features seen in the simulation are very similar to 
those in our observation, both in terms of width and apparent sur¬ 
face brightness excess. We find that regions of high density/surface 
brightness correspond to regions of relatively low magnetic field in 
the map. The thin regions of high density are surrounded by wider 
bands which are highly magnetized. Here, the increased magnetic 
pressure has pushed gas out, decreasing its density. This opens up 
the intriguing possibility that the observed bright bands may be due 
to denser ICM layers between bands of more magnetized plasma, 
where the magnetic pressure is amplified to ~ 5-10 per cent of the 
thermal pressure. Though the magnetic fields in this simulation are 


strong, they do not completely prevent the development of KHI, 
consistent with the indications of KHI in our observation. 

By extrapolating the pressure profiles outside of cold fronts 
inwards, up to the edge, [Reiss & Keshet| ( |2014| ) infer that in the 
inner parts of the spirals most sloshing cold fronts show a pressure 
jump of Pin/Rout ~ 0.8, which may correspond to the amplification 
of magnetic fields mainly below the discontinuity. The Virgo cold 
front does not show any significant pressure jump, which may be 
due to the amplification of the magnetic fields both above and be¬ 
low the front. Fig. shows that such a scenario does occur in our 
simulation. The magnetic fields above the cold front get amplified 
due to the shear in the large-scale velocity field produced by the en¬ 
counter with the subcluster. The same large scale gas motions then 
draw the amplified magnetic fields inward, towards the cold front. 

The observed quasi-linear features may potentially also be due 
to compressed thermal gas formed in front of bubbles entrained 
by the magnetized sloshing gas, spiraling outwards below the cold 
front. Such bubble entrainment was previously suggested by jKeshet] 
( |2012| ). Future numerical simulations, that include both jet infiated 
bubbles and sloshing, will test this possible interpretation. 

Alternatively, for geometries where the sloshing occurs along 
our line-of-sight, KHI would develop predominantly in the plane of 
the sky, forming rolls projected on the bright side of the front. For 
such geometries, Kelvin-Helmholtz rolls projected on the sloshing 
surface could, in principle, appear as quasi-linear features. How¬ 
ever, based on the comparison of the numerical simulations to the 
observations, such geometry is unlikely. None of the simulations 
that we performed show KHI in projection that would look like the 
observed quasi-linear features. 


5.6 The flow of the sloshing gas 

Under the cold front interface, the ICM entropy is lower and above 
the interface it is higher than the radial average. The low entropy 
gas under the interface also has roughly two times higher metallic- 
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ity than the gas above the interface. This entropy and metallicity 
distribution indicates convergent gas flows at the cold front inter¬ 
face, with the low entropy metal rich gas moving outward and, in 
the same time, the high entropy gas moving inward. Such conver¬ 
gent flows are also seen in numerical simulations (e.g. |Ascasibar &| 
|Markevitch|2006[[Roediger et al.|201l| t. 

The entropy profile shows a remarkably flat distribution under¬ 
neath the cold front, extending across a radial range of ~ 50 kpc. 
The flat entropy distribution may indicate that, in this region, the 
gas uplifted by sloshing originated at the same radius. Similar en¬ 
tropy plateaus have also been seen in e.g. the Perseus and Abell 
2142 clusters and interpreted as evidence for gas uplift by slosh¬ 
ing ( [Simionescu et al.|2012[|Rossetti et al.|2013| l. However, if the 
distribution of the uplifted gas deviates significantly from spherical 
symmetry centered on M87 (as assumed by our deprojection analy¬ 
sis), then the shape of the entropy profile measured underneath the 
cold front may be biased. 


6 CONCLUSIONS 

We have analyzed a new, very deep (500 ks) Chandra observation, 
together with complementary archival XMM-Newton data, and per¬ 
formed tailored numerical simulations to study the nearest cluster 
cold front in the sky. We And that: 

• The northern part of the front appears the sharpest, with a 
width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per 
cent confidence). Everywhere along the front, the temperature dis¬ 
continuity is narrower than 4-8 kpc and the metallicity gradient is 
narrower than 6 kpc, indicating that diffusion, conduction and mix¬ 
ing are suppressed across the interface. Such transport processes 
can be naturally suppressed by magnetic fields aligned with the 
cold front. We also And that the temperature jump across the cold 
front is inconsistent with the temperature jumps produced in sim¬ 
ulations with anisotropic thermal conduction, indicating that con¬ 
duction is also suppressed along the magnetic held lines. 

• However, the northwestern part of the cold front is observed to 
have a nonzero width. The broadening is consistent with the pres¬ 
ence of projected KHI eddies on length scales of a few kpc. Based 
on comparison with simulations, the presence of KHI would imply 
that the effective viscosity of the ICM is suppressed by more than 
an order of magnitude with respect to the isotropic Spitzer-like tem¬ 
perature dependent viscosity. 

• Underneath the cold front, we observe linear features that are 
~ 10 per cent brighter than the surrounding gas and are separated by 
~ 15 kpc from each other in projection. Comparison to tailored nu¬ 
merical simulations suggests that the observed phenomena may be 
due to the amplification of magnetic fields by gas sloshing in wide 
layers below the cold front, where the magnetic pressure reaches 
~ 5-10 per cent of the thermal pressure, reducing the gas density 
between the bright features. 

• The deprojected entropy distribution and the metallicity of the 
gas indicate convergent gas flows at the cold front interface, with 
the low entropy metal rich gas moving outward and, in the same 
time, the high entropy gas moving inward. 
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